Dynamics of the particle - hole pair creation in graphene 
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The process of coherent creation of particle - hole excitations by an electric field in graphene 
is quantitatively described. We calculate the evolution of current density, number of pairs and 
energy after switching on the electric field. In particular, it leads to a dynamical visualization of the 

universal finite resistivity without dissipation in pure graphene. We show that the DC conductivity 
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of pure graphene is ? ^-rather than the often cited value of — %-. This value coincides with the AC 

conductivity calculated and measured recently at optical frequencies. The effect of temperature and 

random chemical potential (charge puddles) are considered and explain the recent experiment on 

suspended graphene. A possibility of Bloch oscillations is discussed within the tight binding model. 

PACS numbers: 81.05.Uw 73.20.Mf 73.23. Ad 



c 

o 
o 



> 



o 

On 
O 



X 



1. Introduction. It has been demonstrated recently 
that a graphene sheet, especially one suspended on leads, 
is one of the purest electronic systems. Electronic mo- 
bility reaches values of 2 ■ 10 5 cto 2 / (Vs) and might be 
yet improved [l|, Q indicating that transport in samples 
of submicron length is most likely ballistic. In a simpli- 
fied model of a single graphene sheet (neglecting scat- 
tering processes and electron interactions) the chemical 
potential is located right between the valence and con- 
ductance bands and the Fermi " surface" consists of two 
Dirac points of the Brillouin zone A lot of effort 
has been devoted to the question of transport in pure 
graphene due to the surprising fact that the DC conduc- 
tivity is finite without any dissipation process present. A 
widely accepted value of the "minimal conductivity" at 
zero temperature, 
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was calculated very early on using the Kubo formula in 
a simplified Dirac model as well as in the tight binding 
model 0, H, @, Q . Within this approach one starts with 
the AC conductivity and takes a zero frequency limit 
typically with certain " rcgularizations" (like finite tem- 
perature, disorder strength r\ etc.) made and removed 
at the end of the calculation. As noted by Ziegler Q 
the order of limits makes a difference and several other 
values different from o\ were provided for the same sys- 
tem. The standard value <j\ is obtained using a rather 
unorthodox procedure when the DC limit to — > is made 
before the zero disorder strength limit 77 — > is taken. If 
the order of limits is reversed one obtains [7f 
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When the limit is taken holding lu = rj one can even ob- 
tain a value of 0-3 = ir^ Q, thus solving the "missing 7r" 
problem. Indeed, at least early experiments on graphene 



sheets on Si substrates provided values roughly 3 times 
larger than o\ [||. Recent experiments on suspended 
graphene 0] demonstrated that the DC conductivity is 
lower, 1.7<7i, as temperature is reduced to AK. Hence 
one still faces the question of what is the proper theo- 
retical value. Since the conductivity of clean graphene 
in the infinite sample is a well defined physical quantity 
there cannot be any ambiguity as to its value. 

In contrast both the experimental and the theoreti- 
cal situation for the AC conductivity in the high fre- 
quency limit is quite different. The theoretically pre- 
dicted value in the Dirac model is 02 inde pen dent of 
frequency under condition u) » T/H 0, Tioj ] - The 
Dirac model becomes inapplicable when u> is of order 
of 7/fi, = 4 • 10 15 Hz or larger, where 7 = 2.7eV is the 
hopping energy of graphene. It was shown theoretically 
using the tight binding model and experimentally in [3j 
that the optical conductivity at frequencies higher or of 
order 7//1 becomes slightly larger than 02- Moreover, in 
light transmittancc measurements at frequencies down 
to 2.5 • 10 15 Hz it was found equal to cr 2 within 4%. The 
model does not contain any other time scale capable of 
changing the limiting value of AC conductivity all the 
way to oj— > 0. Therefore one would expect that the 
DC conductivity even at zero temperature is a% rather 
than (j 1 . As we show in this note this is indeed the case. 

The basic physical effect of the electric field is a 
coherent creation of electron - hole pairs mostly near 
Dirac points. To effectively describe this process we de- 
velop a dynamical approach to charge transport in clean 
graphene using the " first quantized" approach to pair cre- 
ation physics similar to that used in relativistic physics 
flil ]. To better visualize the phenomenon of resistivity 
without dissipation, we describe an experimental situ- 
ation as closely as possible by calculating directly the 
time evolution of the electric current after switching on 
an electric field. In this way the use of a rather formal 
Kubo or Landauer formalism is avoided and as a result 
no regularizations are needed. The effects of temperature 
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and of charge fluctuations or "puddles" are investigated 
and explain the temperature dependence of conductivity 
measured recently in suspended graphene Q. Although 
we consider an infinite sample the dynamical approach al- 
lows us to obtain qualitative results for finite samples by 
introducing time cutoffs like ballistic flight time. Various 
other factors determining transport can be conveniently 
characterized by time scales like the relaxation time for 
scattering of phonons or impurities. 

2. Time evolution of the current density at zero tem- 
perature. Electrons in graphene are described sufficiently 
well by the 2D tight binding model of nearest neighbour 
interactions 0], namely with (second quantized) Hamil- 
tonian being sum over all the links connecting sites on 
two triangular sublattices A, B. The Hamiltonian in mo- 
mentum space is 
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-jJ2a elkda w ith 7 being the hop 
ping energy; 8 1 = f (0, V3) , £ 2 , 3 = f(±|,-^ 
are the locations of nearest neighbours on the hon- 
eycomb lattice separated by distance a ~ 3A In 
the Brillouin zone of the lattice there are two Dirac 
points K. = ^(i^), K + = % (|,0) in which 
the energy gap between the valence and the conduc- 
tion band vanishes. Expansion around K-, h{k) = 
Hv g exp (—*f) (Afc x + iAky) , where the graphene veloc- 
ity is v g = leads to relativistic equations for the 
Weyl field constructed as tp 1 = ij\ , -02 = e ~^i ) 2 ■ 

Let us first consider the system in a constant and ho- 
mogeneous electric field along the y direction switched 
on at t = 0. It is described by the minimal substitution 
p = ftk + | A with vector potential (choosing a gauge in 
which the scalar potential is zero) A = (0, — cEt). Since 
the crucial physical effect of the field is a coherent cre- 
ation of electron - hole pairs mostly near Dirac points 
a convenient formalism to describe the pair creation is 
the "first quantized" formulation described in detail in 
[l2T ]. The second quantized state at T = which evolves 
from the zero field state in which all the negative en- 
ergy (— \h (k)\) states are occupied is uniquely character- 
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obeying the matrix Schrocdinger equation ihd t ip = H-ijj 
in sublattice space with the initial condition 



ized by the first quantized amplitude ip k (t) 
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Here Uk is found as a negative energy solution of the time 
independent Schrocdinger equation prior to switching on 
the electric field, Huk = — \h\ Uk. 



A physical quantity is usually conveniently written in 
terms of tp. For example the current density (multiplied 
by factor 2 due to spin) is 
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To first order in electric field ipk = e«'' l ' t (uk + E£ k + ...) 
and consequently J y = Jo + Ea, where 
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The solution of the Schrocdinger equation for the correc- 
tion £ fc is 
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where v k = ^ fft*/|ft|J' Substituting this into Eq. 
the conductivity becomes 
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The zero field current J and the first term (linear in 
time) in the conductivity vanish upon integration over 
the Brillouin zone, since one can choose it to be J gz = 

f n ^ a , dk T f , .„ \ dk v and the integrand is a deriva- 

J-TT/a x J_27r/(3 1 /2 Q ) V b 

tive of a periodic function. 

The integral of the second part (oscillatory in time) of 
Uk (t) gives the result shown in Fig.l. After an initial 
increase over the natural time scale t 1 = ft/7 = 2.5 • 
10~ 16 s, it approaches CT2, Eq.®, via oscillations. The 
amplitude of oscillations decays as a power 
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for t >> t-y. The limiting value is dominated by contri- 
butions from the vicinity of the two Dirac points in the 
integral of Eq.©. The contribution of a Dirac point is 
obtained for t >> t 7 by integrating to infinity (in polar 
coordinates centered at the Dirac point) 
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FIG. 1: Evolution of the current density a (t) = J (t) /E after 
a DC electric field is switched on at t = 0. Unit of time 
is t~, = ft/7. Conductivity is compared to its "dynamical" 
value o 2 = |p The zero temperature conductivity (blue) 
approaches 01, while finite temperature depresses the pair 
creation and eventually the current density vanishes as i . 
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does not influence the result. 

A physical picture of this resistivity without dissipa- 
tion is as follows. The electric field creates electron - 
hole excitations in the vicinity of the Dirac points in 
which excitations are massless relativistic fermions. For 
such particles the absolute value of the velocity is v g and 
cannot be altered by the electric field and is not related 
to the wave vector k. On the other hand, the orienta- 
tion of the velocity is influenced by the applied field. The 
electric current is ev, thus depending on orientation, so 
that its projection on the field direction y is increased 
by the field. The energy of the system (calculated in a 
way similar to the current) is increasing continuously if 
no channel for dissipation is included. Obviously at some 
time the system goes beyond linear response into Bloch 
oscillations which are briefly discussed below. We have 
performed a similar calculation for the evolution of the 
current density for an AC electric field switched on at 
t = 0. After a short transient one obtains the value of 
the DC conductivity 02 independent of frequency. This 
is consistent with both the Kubo formula derivations [Io| 
and optical experiments 

3. The temperature dependence and effect of charge 
"puddles" . At finite temperature T within the first quan- 
tized formalism one adds the contributions of all the ener- 
gies including positive ones weighted with the Boltzmann 
factor. Due to electron - hole symmetry the contribution 
to conductivity of a positive energy state with momen- 
tum k is minus that of the contribution of the negative 
energy state with the same wave vector. This results in 



FIG. 2: The minimal conductivity as a function of tempera- 
ture for time thai = 500t 7 is compared with measured in the 
0.5 jj,m long sample in ref.0]. Values for the random Fermi 
energy are also taken from ref.Q]. 



the thermal factor 
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The first term still vanishes, while the second gives 
a depressed value compared to that at T = 0, see Fig. 
1. Moreover, the conductivity vanishes at the large time 
limit. This is easy to appreciate qualitatively: the contri- 
butions from the vicinity of the Dirac points, \h\ << T . 
which were the main contributors to a (T) are effectively 
suppressed. Physically this suppression can be under- 
stood as follows. As mentioned above the finite resistiv- 
ity of pure graphene is due to pair creation by an electric 
field near Dirac points. The pair creation is maximal 
when in the initial state the valence band is full and the 
conductance band is empty. Thermal fluctuations create 
pairs as well. In the formalism we adopted the finite tem- 
perature initial state is described by the density matrix 
which specified the number of incoherent pairs present in 
the energy range near the Dirac points. Therefore pair 
creation by an electric field is less intensive due to the 
diminished phase space available and the conductivity 
vanishes at large times. 

Under assumption of Dirac point dominance, T << 7 
(definitely covering the temperature range T < 200K be- 
yond which scattering is not negligible [l|), the expression 
can be simplified in the same way as Eq. (fTO]) . 
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and is a monotonically decreasing function of the product 
tT. For t » ty t a T (t) = x W- 

Assuming ballistic transport in a finite sample of sub- 
micron length determining an effective ballistic time tb , 
this contribution cannot explain the increase of conduc- 
tivity with temperature in suspended graphene reported 
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in 0]. However, there is an important source of pos- 
itive contribution to conductivity even in the ballistic 
regime. It was clearly demonstrated that a sample close 
to minimal conductivity consists of positively and nega- 
tively charged puddles. This means effectively that even 
at minimal conductivity the chemical potential fj, locally 
is finite, rather than zero, albeit small on average. Phys- 
ically this implies that in addition to the novel constant 
contribution due to pair creation, there is an ordinary 
contribution due to acceleration of electrons like in or- 
dinary metal. In ballistic regime it grows linearly with 
time. 

The experiment shows that the amplitude of the 
random Fermi energy increases linearly with temperature 
Mt = Mo + a T- For example, for the 0.5/xm long sample 
fi Q = 8meV and a = O.lmeV/ K . The difference between 
and c M= o is equal to the integral in Eq. ([6]) over the two 
regions around the Dirac points determined by \h (k)| < 
fi. That way one obtains for t » t 7 
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which is a monotonically increasing function of the prod- 
uct \xt only (Si is the sine integral function). In Fig. 2 we 
fit the value of the ballistic effective time £& a ; ~ 2- 10~ 13 s, 
which is of the same order of magnitude as for the 0.5/im 
long sample, L/v g ~ 5 • 10~ 13 s. 
4. Discussion and summary. 

To summarize, we studied the dynamics of the parti- 
cle - hole pair creation by calculating the time evolution 
of current density, particle - hole number and energy af- 
ter the electric field is switched on. After a brief tran- 
sient period (of order of several t~ = the current 

density approaches a finite value. The minimal DC clec- 

2 

trie conductivity at zero temperature is different 

from an accepted value The later value was ob- 

tained for nonidcal systems by taking various limits (im- 
purity strength etc.) or in theory of finite size effects [l3| 
and does not characterize an ideal pure infinite graphene 
sheet. At finite temperature T the current density di- 
minishes on the scale of tr = H/T = ^t 1 . Therefore the 
phenomenon of finite resistivity without dissipation dis- 
appears unless there exists a shorter time scale intercept- 
ing the process like 2ir/ui for AC field, relaxation time r 
for scattering off impurities or phonons or ballistic flight 
time tb a i for finite samples. The effect of small random 
chemical potential was also considered. 

Let us now address the issue of the validity of the linear 
response approximation used. Since the model does not 
provide a channel of dissipation, the problem is nontriv- 
ial. Where does the Joule heat aE 2 go? The dynamical 
approach allows us to calculate the evolution of energy 



as well as to go beyond linear response. Of course the en- 
ergy continuously increases with time and at certain time 
approaches the conduction band edge at which stage lin- 
ear response breaks down. We calculated the evolution 
of current density, energy and pair number beyond linear 
response and found that Bloch oscillations set in with a 
period of t B i oc h = = 5^*7- The range of appli- 
cability of the linear response was also determined. The 
average current over larger times is zero. This means that 
at very high fields the minimal conductivity phenomenon 
disappears. However in order to reach the conditions for 
observation of the Bloch oscillations in graphene all other 
time scales r, tT,tb a i, 2tt/lu should be larger than tBioch- 
Additional phenomena beyond linear response as well as 
their relation to the Schwinger's calculation of the pair 
12l , 14j ] is under investigation. 
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